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Abstract : We consider the possibilities of testing CP conservation in the process 
e + e~ — > ZH. The easiest is to measure the forward-backward asymmetry of the Z boson, a 
nonzero value of which signals CP violation. The effects of CP violation are predicted in the 
two Higgs doublet model. We show here that a choise of optimal observables can improve 
the possibility of observing CP violation. Although the number of events is not large, useful 
information about CP violating parameters can be extracted. 



1. Introduction 

After its discovery in the Kaon system 30 years ago[l], CP violation is still a puzzle. 
The standard model can accomodate CP violation in the Kaon system, where CP violation 
is introduced through the imaginary part of the CKM matrix elements. However, because 
of the special structure of the CP violating interaction[2] in the standard model, predicted 
CP violating effects in high energy processes, such as high energy scattering, are extremely 
small. An explicit example is given in [3]. Hence any CP violation found in high energy 
processes will indicate the new physics beyond the standard model. In any case, the 
standard model is generally regarded as a low energy effective theory which will be modified 
or extended at high energies. One of the possible extensions is the two Higgs doublet model 
with two Higgs doublets instead of one. In this model CP violation can arise not only from 
the complex CKM matrix elements, but also from the Higgs sector [4], where, for example, 
the neutral Higgs bosons with different CP eigenvalues can mix. CP violation from this 
mixing may have sizeable effects in high energy scattering processes, which are studied 
in [5,6,7,8]. This model also allows a significant electric dipole moment for elementary 
particles [9]. 

In this work we will study possible CP violation from the two Higgs doublet model 
in the process e + e~ — > ZH, where the initial beams are unpolarized. In this process it 
turns out that a CP test is possible without the knowledge about the polarization of the 
final state. A nonzero forward-backward asymmetry of the Z or H boson indicates CP 
violation. This is in contrast to the case of e + e~ — ► particle + antiparticle, where a CP 
test is only possible when one measures the polarizattion of the final state. However, we 
will also consider the possible information about CP violation revealed by the polarization 
of the Z boson, as the leptonic decay of Z can measure its polarization. 

In this work we also consider the optimization of CP violating effects. In experiment, if 
the number of events available for CP test is not large, the statistical error can be the main 
source of error in the observation. It is possible to reduce the statistical error by a suitable 
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optimization procedure. A general method for optimization was proposed recently [10, 11]. 
We apply the procedure to the present process and compare the effects with and without 
optimization. 

We begin in Sect. 2 with our general formalism and introduce some simple CP odd 
observables. In Sect. 3 we estimate the effect of CP violation from the two Higgs doublet 
model and the sensitivity of our simple CP odd observables to CP violating parameters 
in the model. In Sect. 4 we apply the optimisation procedure mentioned above. Sect. 5 
includes a discussion and our conclusion . 

CP violation in e + e~ — »■ ZH is also discussed in [12] in relation to CP violation in the 
Higgs decay, where only the real part of possible CP violating form factors is considered. 
In this case a nonzero forward-backward asymmetry can not be produced. 

2. CP odd observables 

We consider the following process with the unpolarized initial state: 

e+(p+) + e-(p_) -► Z(k) + H(k H ) -► £ + (k+) + C (fc_) + H(k H ) (2.1) 

The momenta of the particles in the C.M.S. of the initial state are indicated in the brakets. 
The lepton pair comes form the Z boson decay. At the tree level, the Higgs boson in the 
process (2.1) is produced through Bremsstrahlung off the virtuel Z boson produced by 
e + e~ annilation in the standard model. 

Instead of the momentum fc+ in Eq. (2.1), we will use the momentum q + of £ + in the 
Z rest frame, which is of course related to k_|_ by a Lorentz boost. We use the notation k to 
indicate a unit vector in the direction of k. We can introduce a density matrix R(p+, k, q + ) 
for the process in Eq. (2.1), which we normalize so that 

a(e + e~ ^ ZH ^e+e~H) 1 |k| 1 f 1 



J ( m-±-Jdn + R( P+ ,k 1 c l+ ) (2.2) 



Br(Z^e+£-) 2s 47rv/i 4tt 

Here s = (p+ +p_) 2 and dQ(dQ+) is the solid angle of the vector k(q+). The quantity R 
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can always be written in the form: 

#(P+,k,q + ) = A(p + ,k) + q + -B(p+,k) + (q +i q +j - -S ij )C ij (p+,'k) (2.3) 

If CP invariance holds, we have: 

A(p + ,k) = A(p+,-k), B(p+,k) = B( P+ ,-k), C ij (p+,k) = C ij (p+,-k) (2.4) 

To test CP conservation one therefore has to verify the equations in Eq. (2.4). In Eq. (2.3) 
the quantities B and Cy contain the information about the polarization of the Z boson. 
Without measuring the Z polarization one can still test CP invariance by checking the CP 
constaint for A in Eq. (2.4). We propose the following CP odd observables for a CP test: 

1 =x = p+k, 2 = (q+-p+)(q+-k), 3 = (q + • p + )q + • (p + x k) (2.5) 

A nonzero expectation value of any observable Oj in Eq.(2.5) signals CP violation. We 
will call these observables simple CP odd observables. Note that the < 0\ > can only 
be nonzero if the quantity A does not satisfy its CP constraint. The forward-backward 
asymmetry of the Z boson can be defined through 0±: 

where N(0\ > 0) and N(0\ < 0) is the number of events number with 0\ > and 0\ < 
respectively. 



3. CP violation from the two Higgs doublet model 

We consider here a two Higgs doublet model as a possible extension of the standard 
model. In this model CP can be violated due to the complex expectation values of the Higgs 
doublets. However, in general CP violation occures together with the presence of flavour 
changing neutral currents at the tree level. The flavour changing neutral currents can be 
eliminated by imposing some discrete symmetry on the model. This dicrete symmetry is 



softly broken in the Higgs potential[4]. In this way it is possible to have CP violation in 
the model with the absence of the flavour changing neutral currents at the tree level. 

In this model, there are three neutral Higgs bosons, two of them have the CP eigen- 
value +1, and the other one has —1 as its CP eigenvalue, if CP invariance holds. In this 
case, the neutral Higgs boson with the CP eigenvalue —1 can not be produced in the pro- 
cess (2.1) at the tree level, because it decouples from two Z bosons. If CP is violated, the 
neutral Higgs bosons do not have any definite CP eigenvalue, but are mixtures of the CP 
eigenstates. In the process we consider the effect of CP violation can occur at one loop 
level only through the type of the diagrams given in Fig.l, where the loop is a fermion 
loop. If the mass of the light flavours is neglected compared to mt, only the top quark 
contributes. In the two Higgs doublet model we write the relevant interaction between the 
neutral Higgs boson H and the top quark and the interaction between two Z boson and 
the Higgs boson H as: 



eM z 
2 sin 6w cos 6 Z 



L H u = -iB t i l5 tH, L ZZH = n „ ; _„ 1„„ d ■ Z^Z^H (3.1) 



The parameters Bt and d are real. Furthermore \d\ < 1. If Bt and d are simultaneously 
nonzero, CP is violated. 

We parameterized the CP violating part of the VHZ vertex with a form factor Fy as 

z sin u\y cos u\y 

where the indices \x and v are the Lorentz indices for the polarization of the virtual vector 
boson V and the real Z boson respectively. V stands for a virtual photon or a virtual 
Z boson and has the momentum q with q 2 = s, where s > [Mh + Mz) 2 - The on-shell 
Z boson carries the momentum k. Calculating the relvant diagrams, one finds that the 
real part of the form factor Fy has a complicated form, expressed with Spence functions. 
Since we will take y/s = yep = 500GeV, and the mass of the top quark will be taken as 
175GeV, as suggested by recently in experiment [13], the absorbtive part of the vertex, i.e. 
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the imaginary part of the form factor, is always nonzero in the s channel. We calculate 
the imaginary part first, the real part can then be obtained from the dispersion relation. 
For V = Z the imaginary part is: 

sin t)\y COS t>w 

fy <«> = IwXyT* ' ^W = ^^«^-^ + /«^)^ + 2A^ (3.3) 
_ (s - Ml - M\) -( S -M 2 H + M 2 z )(3 t p z 
L {s-M 2 H - M 2 Z ) + ( S -M 2 H + M z )(3 t (3 z 

Here fit and (3z are the velocity of the top quark and the Z boson in the rest frame of the 

virtual Z, and Ez is the energy of the real Z boson. Im_F 7 (s) can be obtained from Eq. 

(3.3) by the substitution: (g v ) 2 — > | sin 6w cos 9w9v an< ^ (9a) 2 ~^ 0- From the imaginary 

part we can obtain the real part: 

ReFy(s) = -P / cfe' Vl J , so = max{4m 2 (M H + M z ) 2 } (3.4) 

Note that there is also another contribution to the imaginary part of the form factor, if 
the mass of the Higgs boson is larger than 2mt. We will not consider this case. The 
CP violating part of the density matrix R is generated by the interference between the 
amplitude at tree level and the amplitude including the vertex in Eq.(3.2). The result for 
Ris: 

ed \ 3s 



* (p+ ' k ' q+)= U a W»J W^MW m{l ^ 2 ) + ^+-^ 



+ 2Mz(E z - M z )x(p+ • q+)(k • q+) + (k • q_ + Y(-(3 z E z + x z (E z - Mz) 1 } 

+ ^j^Pz sm9 w cose w (l - ^)hmF 1 (s)[(p+ • q+)(k • q+)M z E z 
Mz s a 

+ (k • ci + ) 2 xMz(M z - E z ) + xM\\ 

- y^ReF z (s)PzE z q+ • (p+ X k)[M z ( P+ • q + ) + x(E z - M z )(k • q + )]| 

(3.5) 
Here we used the fact that the coupling constant g v (— — \ + 2 sin 2 9w) of the Z boson to 
the neutral vector current is much smaller than the coupling constant g A to the neutral 
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axial current and neglected g v . The CP conserving part of R at the tree level is also given 
in Eq. (3.5). With the density matrix we can now calculate the expectation value of the 
simple CP odd observables and the CP asymmetry in Eq. (2.6). These quantities are 
proportional to ^-: 

<O t >=Wi-^, a C p = y—t (3.6) 

Further, < 0± t 2 > receive contributions only from the imaginary part of the form factors, 
while the real part contributes to < O3 >. We take y/s = 500GeV and mt = 175GeV and 
obtain the numerical results for the coefficients in Eq. (3.6): 

wx = 5.6 x 1(T 3 , w 2 = 3.8 x 10" 3 , w 3 = -4.4 x 1(T 4 , for M H = lOOGeV 

(3.7) 
wi = 6.7 x 10" 3 , w 2 = 4.3 x 10" 3 , w 3 = -4.6 x 10" 4 , for M H = 200GeV 

For the density matrix R given in Eq. (3.5), y has a simple relation to w±: y = |wi- The 
result for j/asa function of the Higgs mass Mn is plotted in Fig. 2. From Fig. 2 note that y 
increases as Mh increases, mainly due to the decreasing cross section with the increasing 
Higgs mass. 

To study how sensitive our observables are to the CP odd parameter ^ , we introduce 
the sensitivity pi of a observable Oi : 



P. = 44^ (3.8) 

\wA 



Remembering that the statistical error of a observable Oi is given by: 



where N is the number of the available events for measuring the observable Oi, the sensi- 
tivity pi has the meaning that the parameter -j- must be at least larger than -4= to have 
a measurable effect on Oi. The results for pi are: 

Pl = 88, p 2 = 80, p 3 = 558 for M H = lOOGeV 

(3.10) 
Pl = 75, p 2 = 72, p 3 = 532 for M H = 200GeV 
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Pi as a function of Mh is again plotted in Fig. 3. From Eq.(3.10) we see that the observable 
01,2 are much more sensitive than O3. The sensitivity for a>cp is given simply by y _1 , and 
it is about 30% larger than p\. 

4. Optimization 

In [9,10], it is discussed how to construct the most sensitive observable to the value of a 
physical parameter in a reaction, where the distribution is linear in the physical parameter. 
We give here briefly the result from [10]. Let us consider a scattering process where the 
differential cross section has the form: 

E(0)#=(E o (0)+0-£i(0))d0 (4.1) 

Here <p denotes the sample of the relevant phase-space variables, with which we construct 
observables. The observable T which is most sensitive to the parameter g in Eq. (4.1) 
takes form: 

T = T ^ = ^§ ^ 

This result is generalized in [11] for more complicated cases. Such types of the observables 
may be called optimal opservables, following [11]. 

In our case, with the distribution, i.e. the density matrix R given in Eq.(3.5) we 
construct according to Eq. (4.2) the optimal observables for the form factors in Eq.(3.2) 
which we want to measure. The observables are: 

(p+ • q+)(q+ • k)M z E z + (q+ • k) 2 xM z (M z - E z ) + xM 2 
±0 — 

Do 

_ -f3 z E z q+ ■ ( P+ x k)(M z (p + ■ q + ) + x(E z - M z )(q + ■ k)) 
Ts ~ W (4-3) 

D Q = E\{\ - (3 2 z x 2 ) + M 2 (p+ ■ q+) 2 + 2M Z (E Z - M z )x(f>+ ■ q+)(q + ■ k) 
+ (q+ • k) 2 (-(3 2 z E 2 z + x 2 (E z - M z ) 2 
The optimal observables T<i and T3 correspond to the simple observables O2 and O3. To 



present the numerical results we define in anaolgy to Wi and pi the quantities: 



<r,>=3, pt, = ^^ (4.4) 

CI \Zi\ 

The numerical results are: 

2 2 =4.2xl0" 3 , z 3 = 1.0 x 10" 3 , Pt 2 =67, p Tz = 401 for M H = lOOGeV 

(4.5) 
z 2 = 4.2 x 10" 3 , z 3 = 9.3 x 10" 4 , p n = 60, p n = 414 for M H = 200GeV 

Comparing the sensitivity given in Eq.(3.10) for the simple observables the sensitivity of 
the optimal observables is from 20% to 40% better. However, the optimal observable T 3 
is much less sensitive than T 2 , and even much less sensitive than Oi ; 2- 

To construct the optimal observable T\ corresponding to Oi, we integrate out the 
solid angle of the letpton momentum, i.e. extract A(p+,k) in Eq. (2.3): 



p 2 rt 2 F 2 q 



^(P+, k) = (—5 — Y- , z 9NO • <^ 2 - /3i - /3ix 

VH+ ' ; sin 2 ^ w cos 2 0vj/ 32(s-M|) 2 \ Mz Mz 

+ 8(1 - ^) sin 6 W cos e w ^^^- lmFJs) ■ x 
s E z 



(4.6) 



In this simple case we find that the effect of the improvement with the optimal observable 
is tiny, as also claimed in [10]. Hence we do not present our numerical results for 7\. 

5. Conclusion 

In this work we studied CP violating effects from the two Higgs doublet model in 
e + e~ — ► ZH. The effects depend on a unknown parameter -4- of the model, defined 
in Eq. (3.1). To detect these effects we propose two sets of CP odd observables. One 
set contains the simple CP odd observables, the other one contains the so-called optimal 
observables, which are designed specially to maximize the CP violating effects from the two 
Higgs doublet model, or more generally from the form factors in Eq. (3.2). The sensitivity 
of these observables to the parameter -4- is estimated. From our results, however, the 
optimal observables do not improve the sensitivity to the CP violating effects very much, 



the improvement is about 20% to 40%, the price for this improvement is that the optimal 
observables take a more complicated form than the simple observables. 

One of the simple observables, i.e. the forward-backward asymmetry acp of the Z 
boson or the corresponding mean value < 0\ > should be readily accessible in in experi- 
ments. In principle, since the lepton momentum is not involved here, more events should 
be available for measuring 0\ and a>cp than for the other observables 02,3 an d ^2,3. Com- 
paring the results obtained in the previous sections, we conclude that the observable Oi 
or acp are at most sensitive to the CP odd parameter -4-. However, the other observ- 
ables should also be studied in experiments, because our conclusion is specific to two Higgs 
doublet model, and it is possible in other possible extensions of the standard model that 
< Oi > and acp both vanish, but others do not vanish. 

Finally, we note that one can determine how large Bt should at least be, independent 
of the parameter d, to have a possiblly observable effect, if one knows the luminosity of the 
e + e~ collider. Assuming the luminosity to be 10fb~ per year and requiring CP violation 
to be such that the expectation value of 0\ is larger than twice of its statistical error, then 
we require for Mh = lOOGeV: 

B t > 7.2 (5.1) 

for the observablity of the effect. 
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Figure Caption 

Fig.l: One of the Feyman diagrams contributing to the form factor in Eq. (3.2) 
Fig. 2: The parameter y for the CP asymmetry acp as a function of the Higgs mass. 
Fig. 3: The sensitivity p\ for the observable 0\ as a function of the Higgs mass. 
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